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Abstract 

We examine the nature of the unknown enhancement around 3 GeV observed by the BABAR 
cohaboration in the rupp spectrum of the — )■ ppD^ decay. After demonstrating that the peak 
can be neither identified as a charmonium state, such as r/c or J/ip, nor classified as one of the hght- 
flavor mesons, we conclude that it corresponds to a glueball fitted as X(3020) with {mx, ^x) = 
(3020 ± 8, 103 ± 3) MeV, which could be the first glueball state above 3 GeV. 
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Introduction — The glueball (G) is a bound state that contains no valence quarks but 
gluons only. This is because gluons, which are charged with colors in QCD and force carriers 
to bind quarks becoming mesons and baryons, can also glue themselves together to form a 
bound state. Since it is a unique feature purely for the non-Abelian gauge fields, whether 
the existence of the gluon condensates can be well established or not appears to be a real 
test for QCD. In principle, the search for glueballs depends on gluon-rich processes, such as 
the radiative J/\E' decays via cc — )■ 'ygg. 

However, it is known that t^ 



predicted mass around 1.7 GeV 



le glueball identifications are complicated [liJal. With the 
^,|5|, the lightest J^*" = 0"^"'" scalar glueball is allowed to mix 



with nearby qq mesons in the spectrum. Since there are two states, /oQSOO) and /o(1710), 
proposed to be composed of the glueball in different mixing scenarios , the identification 
is obscure. The lightest tensor glueball with J^*" = 2"*""'" is believed to have a mass close 
to 1.3 GeV in the MIT bag model I?] and 2.4 GeV in the lattice QCD calculation Q, 5|. If 
M(2++) = 1.3 GeV, both /2(1270) and /2(_1525) as the ground states of the 2++ mesons are 



le other hand, if M(2 



2.4 GeV I2I 



argued to have the 2"'"+ glueball content [8|. On t 

both /j(2220) (J = 2 or 4) Q and /2(2340) llll are considered to be the candidates, 
whereas the existence of /j(2220) is not conclusive [12|. Unlike O"*"^ and 2+"*", the difficulty to 
establish the lightest ^ pseudoscalar glueball is that the predicted M(0 ^) ~ 2.6 GeV in 
the lattice QCD calculation 0, 5| has no correspondence with the data. Nonetheless, ?7(1405) 



seems to be a perfect candidate [13[. Particularly, the unseen in 77 reactions [IJ] reflects 
that its components are gluons. In addition, Xfl835) first measured in the J/\E' — )■ 7pp 
decays 15| is another possible glueball state 16| at a mass below 2 GeV. Interestin gly, 



instead of taking the candidates as the pure glueballs, the rj — rj' — G [17[ and rjc — G 
mixing scenarios for ?7(1405) and X(1835) are able to allow their own glueball components 
to be at least 2 GeV, respectively. This alleviates the problem that the glueballs with masses 
smaller than 2.6 GeV. However, due to the two mixing scenarios, it is clearly difficult to 
draw a clear conclusion about the glueball state. Moreover, as the mesons lying between 
the glueball candidates are not necessarily well settled, it increases the uncertainty. As a 
result, it is unrealistic to say that the search for the lowest glueballs can provide a compelling 
evidence for their existences. 

Since the searches for the three lightest glueball states are still unclear, we should move 
on to explore the heavier states. Presently, as the PANDA experiment built to scan heavy 



glueballs with masses under 5.4 GeV will not be ready until 2018, we can only use the decays 
of the charmonium states, such as ?7c, J /ip and ip{2S), in the mass range of 3.0 — 3.7 GeV, 
where the glueballs with masses around 3 GeV have been richly predicted. 

On the other hand, although the B meson decays are not regarded as the gluon-rich 
processes, they can be more beneficial to offer accesses to a wider detection range of heavy 
glueball productions. We note that the three-body baryonic B — )■ ppM decay with a two- 
step process B ^ {G ^ pp)M could be an ideal channel, where M is the recoiled meson. 
In particular, one can think of the G ^ pp transition as an inverse process of the pp 
annihilation, which has been used at LEAR and PANDA as a gluon-rich process to search 
for glueballs. In fact, the process of B ^ — ppK has been applied to constrain the 



narrow resonant state known as the glueball candidate /j(2220) 19|, |20|. Recently, the 



BABAR collaboration has o 



Dserved an unknown enhancement at 3.0 — 3.1 GeV in the nipp 



spectrum of — ?■ ppD^ 2l|. We shall demonstrate that the peak is a sign for a resonant 
state as it is unable to be reproduced by the perturbative QCD (pQCD) calculations. Since 
the charmonium states, such as rj^ and J / tp as well as the light-flavor mesons are not favored, 
we have to introduce the glueball state at a mass above 3 GeV as the resonant state. 

Data Analysis — For the rripp spectrum in the B^ — ppD^ decay, the data around 3 



GeV measured by the BABAR collaboration [2l| are robust as the highest data point is 
about 5a significance away from zero. We emphasize that the peak at 3 GeV can only be 
revealed when the rripp spectrum is split into two parts in terms of rriDp > 3 and < 3 GeV, 
respectively |2l|. We note that the threshold effect that shadows other peaking signals has 
been removed to the rripp spectrum according to rriup > 3 GeV. 

To find out if the peak at 3 GeV can be accounted in the theory without introducing a 



new state, we start with the amplitude based on pQCD counting rules, given by 22 1 



AiB'^ppD'') = ^V,,V:MD\-cu)v-A\0){pp\{db)y_A\B') , (1) 

where Gi;' is for the Fermi constant, Vcb and Vud represent the CKM matrix elements for the 
b — 7- cud transition at the quark level, and {qiq2)v{A) stands for gi7^(75)q'2- For the D meson 
production, we have 

{D'\{cu)y^A\0) = rfDP'', (2) 
where fo is the decay constant of D. The matrix elements for the B^ pp transition are 



mop<3 GeV 
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FIG. 1. Invariant mass spectrum as the function of the invariant mass nipp in — t- ppD^, where 



the data points are taken from Ref. 



2l| , the sohd hne includes the contributions from the resonance 



and pQCD counting rules, and the dashed Hne corresponds to that without any resonance state. 



parameterized as 23 1 



(3) 



where p = Pb — Pp — Pp and q = Pp + Pp with pi {i = B,p,p) representing the momenta of 
the particles, and the form factors fj{(u) = Df.^g_.)/t^ (j = 1, 2, ■ ■ ■ , 5) with t = are in 



accordance with pQCD counting rules 



24| , such that the form of 1 /t" with = 3 is to count 



the number of the hard gluons for the B ^ pp transition 25|. The theoretical inputs for 



'^2 5 /l>5 and Df.tg.) can be refereed to those in Ref. 



22l |. To integrate over m/jp in the range 



< 3 GeV, we use the equation in Ref. [26| for the phase space of the three-body decay. Our 
result is shown in Fig. [H 

As seen in Fig. [T], the dashed line in the rupp spectrum fits with the flatness of the 



non-peaking data points given in Ref. |2l|, which illustrates the suppression featured by 
fj{gj) oc in pQCD. For the data points unlinked by the dashed line, we suspect that 
they are not related to any resonances as the significances are less than 2.0. S.ce the 
fit with only pQCD counting rules fails to reproduce the enhancement around 3 GeV, a 
resonant state at the peak is suggested. As B^ — )■ (M(cc) — pp)D^ is allowed to take place, 
with the mass of M around 3 GeV, J/ip or t]c can be the candidate for the resonance. In 
Eq. (|3]), the B^ pp transition is via B^{bd) — )■ {dd — pp). In pQCD counting rules, one 
needs three hard gluons for the transition: one hard gluon is to speed up d, while the other 



two attach to the valence quarks inside pp. Without being directly related to pp by the hard 
gluons, the dd pair can be bounded as the light-flavor meson M{dd). It is also possible for 
the dd annihilation, such that the multi-gluons are generated to form the glueball G at a 
mass around 3 GeV. Therefore, we have three possibilities: the charmonium M(cc) such as 
J/ip and ?7c, the light-flavor meson M{dd), and the glueball G. 

As a ly-boson exchange process, — (M(cc) — )■ pp)D^ is naturally suppressed. With 
the relation, given by 

^ ( J/^ ^ PP)D^) ~ ^ J/^D^)B{J/^ pp) , (4) 

we flrst test the resonant J/ip state by taking the measured data in Fig. [T]and B{J /ip — )■ pp) 



m 



26| as the inputs. It turns out that B{B^ — )■ J/j/jD^) ~ 10 ^, which strongly disagrees 



with the predicted B{B^ — ^ J/ipD^) ^ 10 ^ 28|, |29| as well as the experimental upper 



bound B{B^ J/ipD^) < 10"^ |26|. In addition, it is stated in Ref. |2l| that the decay 
width T{J /ip) = 93 keV is not consistent with the broad 100-200 MeV in the rripp spectrum. 
Similarly, we also obtain B{B^ — )■ rjcD^) ~ 10~^, which is much larger than the predicted 



B{B^ — )■ TjcD^) ^ 10"^ 29|. Clearly, the resonance cannot be the charmonium. 

For B^ — 7- {X — pp)D^ with X to be M{dd) or G, the relevant amplitude is the same as 
in Eq. ([1]), while the matrix elements of the B^ — > pp transition are given by 

{pp\{db)y.Am = {pp\X)- {X\{db)v.MB') , (5) 

where mx and Tx are the mass and the decay width, respectively. Consequently, the relevant 
amplitude of B^ — j- [X — > pp)D^ now reads 

Ar{B' ^ (X ^ pp)D') = ^K6K>27T T^n{a + h^,)v , (6) 

V2 [t-mj^)+imxTx 

with the constants a and h. We note that, no matter what spin the X particle has, the 
parameterization for the B^ — {X — )■ pp) transition can be factored into a and h. Although 
a and h are in principle energy-dependent, their values can only be slightly changed with 
the deviation for the decay width around 100-200 MeV compared to the energy range at 3 
GeV. Since the parity determination for the X particle is uncertain, we set \a\ = \b\. By 
taking the 6 neighboring data points for the enhancement around 3 GeV as our inputs to 
the combined amplitude A = A{B^ ppD^) + Ar{B° {X ^ Pp)D^), we flt \a\ = \b\, 



the mass, and decay width for the X particle to be 

\a\ = \b\ =4.7 ±0.9, 

(mx, Tx) = (3020 ± 8, 103 ± 3) MeV , 



(7) 



respectively, where the small errors correspond to the accuracies. Our result with the above 
resonance is presented as the solid line in Fig. [1] From the figure, we observe that it can 
fully explain the peak. 

Due to its mass, X(3020) is not likely to be M{dd). In fact, there is no observation of any 
light-fiavor meson heavier than /e(2510) in the literature 26|, and the predicted spectrum 



of the excited mesons does not span above 2.8 GeV [30|]. This agrees with the study of the 
hadronic Regge trajectories 3l|, where the mass limits are given to be (2.86 ± 0.11) and 
(3.10±0.11) GeV for nn and ss mesons, respectively. Moreover, the heavier meson with the 
quark pair inside in the higher state has more decay channels, resulting in a broader decay 
width. Since /6(2510) has its decay width of (283 ± 40) MeV, it is also impossib 



32 



e for the 



33|, the 



heavier M{dd) to shrink the width back to (103 ± 3) MeV. As stated in Refs. 
glueball can be ideally observed in the mass range above 3 GeV, where the productions of 
the light-fiavor mesons are not able to take place. As a result, it is reasonable to recognize 



X(3020) as the glueball. Furthermore, it is promising that X(3020) can be one of the 
glueballs predicted from various QCD models (4l. Isl. I33l-l36l| in Table [U where the 2~+ glueball 



TABLE I. Predicted glueballs around 3 GeV in Refs 
MeV. 



. Q,y, 3-36 



where the units of masses is in 



JP<^ = 2-+ 


1~ 


1+- 


3100 ± 30 ± 150 [4] 
3040 ± 40 ± 150 [5] 
2950 ± 150 [34] 


3200 ± 200 [3^ 
3240 ± 330 ± 150 [35] 
3020 [36] 


2940 ± 30 ± 140 [4] 
2980 ± 30 ± 140 [5] 
3270 ± 340 [M] 



contains 2 gluons, while the 1 and 1"^" ones are allowed to have 3 constituent gluons. Since 
J/ip{l ) mainly decays into ggg, the O — J/tp admixture with O denoting the 1 glueball 
is proposed to provide the solution to the so-called pvr puzzle 37|. Recently, the experimental 
data from the charmonium decays at BES and CLEOc turn out to disfavor this solution 



Nonetheless, one of the original mixing scheme leads to \mQ — mj/^\ < 80 MeV and Tq < 
120 MeV j39|, agreeing with the fits in Eq. ([7]). Finally, it is interesting to point out that 
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the same resonance also appears in — )■ ppD*^ [2lf]. The combination of the two sets of 
data should be statistically more convincing. 

Discussions and Conclusions — The baryonium is an alternative to explain X(3020), 
while the baryon pair BB in the transition are bounded before the pp production. For 



example, it can be the pp bound state as proposed to explain X(1835) 40| except that the 
binding is in a higher state for a heavier mass. The excited N*N* bound state with A^* being 
one of the states A^(1440), A^(1520) and A^(1535) is also possible, provided that it is allowed 
to release energy to turn itself to be pp. As one of the possible resonant baryoniums, the 
A(1520)A(1520) bound state is nearby X(3020), but with more suppression by the ss — uu 
transition compared to the A^(1520)iV(1520) bound state. Due to the lack of the theoretical 
understanding, the discrimination among the baryoniums is not an easy task. Hence, from 
the theoretical consideration, the glueball seems to be the most promising candidate to 
understand the resonance in Eq. ([7]). 

In sum, we have identified the existence of the glueball state at 3.02 GeV based on 
the peak in the rupp spectrum of B^ — )■ ppD for m£,p < 3 GeV observed by the BABAR 
collaboration, which could be the first glueball state above 3 GeV. Explicitly, it has been 
fitted to be A:(3020) with (m^, Tx) = (3020 ± 8, 103 ± 3) MeV. 
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